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Black hole driven outflows have been observed in some dwarf galaxies with active galactic nuclei [24], and likely 
play a role in heating and expelling gas (thereby suppressing star formation), as they do in larger galaxies [11]. 
The extent to which black hole outflows can trigger star formation in dwarf galaxies is unclear, because work in 
this area has hitherto focused on massive galaxies and the observational evidence is scarce [12} {23} [13]. Henize 
2-10 is a dwarf starburst galaxy previously reported to have a central massive black hole[38}|36) [37/39] though 
that interpretation has been disputed since some aspects of the observational evidence are also consistent with 
a supernova remnant [16,9]. At a distance of ~9 Mpc, it presents an opportunity to resolve the central region 
and determine if there is evidence for a black hole outflow impacting star formation. Here we report optical 
observations of Henize 2-10 with a linear resolution of a few parsecs. We find a ~150 pc long ionized filament 
connecting the region of the black hole with a site of recent star formation. Spectroscopy reveals a sinusoid-like 
position-velocity structure that is well described by a simple precessing bipolar outflow. We conclude that this 
black hole outflow triggered the star formation. 


Radio observations of Henize 2-10 using very long baseline interferometry reveal a nuclear, compact, non-thermal 
source with a luminosity of Lp ~ 4 x 10°° erg s~! and a physical size < 3 pc x 1 pc [B6]. High-resolution X-ray 
observations unveil a point source with Ly ~ 10°° erg s~! that is spatially coincident with the compact nuclear radio 
source [37]. There are two possible explanations for these radio and X-ray observations alone - a highly sub-Eddington 
massive black hole (i.e., a low-luminosity AGN) or a very young supernova remnant[36] [16]. However, there are other 
observational results to consider regarding the origin of the nuclear radio/X-ray source in Henize 2-10. We summarize 
these results in Figure|4](Extended Data Table 1) and demonstrate that a highly sub-Eddington massive black hole is 
consistent with all the available observations including new results presented here, while a supernova remnant is not. 


We observed Henize 2-10 at optical wavelengths using the Space Telescope Imaging Spectrograph (STIS) on the 
Hubble Space Telescope (HST). We obtained observations of the central regions of Henize 2-10 with the 0.2"-slit in 
two orientations. The first, referred to as the EW orientation, is centered on the nuclear radio/X-ray source and aligned 
with the filamentary ionized structure between the two bright, extended regions of ionized gas previously identified 
in narrowband Ha imaging with HST (Figure[I). The second, referred to as the NS orientation, is centered on 
the nuclear source and rotated 90° with respect to the EW observation. We obtained high-dispersion observations 
(velocity resolution of ~ 40 km s~*) at both slit positions using the G750M and G430M gratings, which cover the 
strong emission lines of interest (e.g., Ha, [NI], [SII], [OI] and HØ, [OIM], respectively). 


The kinematics of the ionized gas provide evidence for an outflow originating from a nuclear massive black hole. 
First, we detect substantially broadened emission lines at the location of the central source in both slit orientations. 
In particular, the [OI]6300 emission line has a broad component with a full width at half maximum (FWHM) of 497 
km s7! as measured in the EW slit (and 445 km s~! in the NS slit) using a 3-pixel extraction region in the spatial 
dimension. While the [OI] line is too weak to be detected all along the EW filament (or NS slit), it is detected at the 
location of the bright star-forming region ~70 pc to the east with a FWHM = 103 km s~!, much less than that of the 
central source (e.g., see top left panel of Figure[2). The [OII]5007 line is much stronger than [OI]6300 and detectable 
all along the EW filament. The FWHM of the broad component of [OIN] at the location of the nuclear source is 271 km 
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Figure 1: HST optical image of the dwarf starburst galaxy Henize 2-10 (Image credit: NASA/STScI). The inset 
shows a narrowband Ha-+continuum image of the central 6” x 4” region. At the distance of Henize 2-10 (~9 
Mpc), 1” corresponds to ~ 44 pc. White contours indicate compact radio emission detected with very long baseline 
interferometry and mark the location of the central massive black hole [36], which is also detected in X-rays [37]. Our 
HST spectroscopy demonstrates that the black hole is driving an outflow that is triggering the formation of young 
massive star clusters. The main image is 25” (~1.1 kpc) across. 


s~ 1, which is somewhat broader than the median value per spatial pixel along the EW slit (175 km s~ with a standard 


deviation of 53 km s~'). We do not observe such broad emission beyond the location of the central source in the NS slit 
(see top right panel of Figure[2). We emphasize that the line widths at the location of the nuclear source are consistent 
with a low-velocity outflow from a massive black hole, but would be anomalously low for a very young supernova 
remnant with typical FWHMs of thousands of km s~! [26] [5]. 


In addition to broadened emission lines at the location of the central source, we also find Doppler shifted velocities 
along the EW slit that exhibit a coherent sinusoid-like pattern that is relatively smooth, especially in comparison to the 
position-velocity diagram along the NS slit orientation that shows no evidence of a coherent pattern in the Doppler 
velocities of strong emission lines (see Figure [2|and Methods). Moreover, a simple model of a precessing bipolar 
outflow broadly reproduces the observed sinusoid-like velocity pattern along the EW slit that is aligned with the ionized 
filament seen in narrowband Ha imaging with HST (see Methods for a description of the model). Precessing jets have 
been observed in many AGNs, although they are typically found in more luminous quasars and radio galaxies [14} {10}. 
Theories for the origin of precessing jets/outflows include accretion disk warping, jet instabilities, and the presence of 
massive black hole binaries [34] [31]. On the other hand, the coherent velocity pattern we observe over ~150 pc along 
the ionized filament centered on the central source is incompatible with a supernova remnant origin since supernova 
remnants do not drive quasi-linear outflows on such large scales. 


There is also evidence that the black hole outflow is triggering the formation of star clusters in the central region of 
Henize 2-10. HST imaging shows that the ionized filament extends eastward from the massive black hole to a bright 
knot of ionized gas and site of recent star formation located 1.5" (~70 pc) away from the black hole (Figure 3). Given 
that our HST spectroscopy along this filament exhibits a continuous velocity pattern, which can be tracked from the 
black hole to the eastern star-forming region and is well described by a precessing bipolar outflow model, this strongly 
suggests that the outflow driven by the black hole is causally connected to the region of recent star formation. There is 
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Figure 2: Optical spectra and ionized gas kinematics for the central region of Henize 2-10. Top panels show 2D 
continuum-subtracted HST STIS spectra of key emission line regions along the EW and NS slit orientations. The 
middle panels show position-velocity diagrams along each slit position where Doppler shifts are relative to the systemic 
velocity of the galaxy (873 km s~') [21]. The velocity pattern in the EW slit is much more coherent than that in the 
NS slit. The bottom panels show extracted 1D spectra at the location of the nuclear massive black hole. The region 
surrounding the nucleus (0 pc in both slit positions) shows strong broadened emission lines, including [OI]6300 and 
[OMT]5007. 
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also a secondary, blue-shifted peak (offset by 154 km s71) detected in the emission lines at the location of the bright 
star-forming knot, suggesting the outflow is pushing the line-emitting gas clouds and influencing their kinematics 
(Figure B}. The double peaked lines would naturally arise as the outflow intercepted dense gas and primarily pushed it 
in the lateral direction rather than ahead of the flow. The handful of young star clusters associated with the star-forming 
knot have ages of ~4 Myr (see Methods) and are predominately aligned in the north-south direction, consistent with a 
scenario in which they formed from gas moving in opposite directions due to the impact of the black hole outflow [2O]. 
There is also a local peak in the gas density at the location of the star-forming knot with ne > 104 cm~? (see Figure 
uy consistent with the outflow compressing gas clouds and enhancing star formation in this region. It is notable that 
we also detect double peaked emission line profiles to the west of the massive black hole at the boundary of a dark 
cloud of gas and dust (region 5, see Figure[8), suggesting the other side of the bipolar outflow is also intercepting dense 
clouds and pushing them perpendicular to the outflow direction. Just to the west of this region (region 6, Figure[8), the 
equivalent width of the Ha emission line indicates stellar ages < 3 Myr, indicating a very early stage of star formation 
in which the infant star clusters have not had enough time to destroy or disperse the dense molecular gas from which 
they formed, leading to high levels of extinction [4]. 


The ionization conditions of the gas provide additional support for a massive black hole with an outflow. First, our 
HST spectroscopy at ~0.1” resolution clearly reveals non-stellar ionization at the location of the central source based 
on the flux ratio of [OI]/Ha, which is consistent with gas photoionized by an accreting massive black hole (Figure 
p). Additionally, supernova remnants are known to be strong [SII] emitters with log([SU]/Ha) > -0.521, which is 
not seen in the STIS spectrum of the nuclear source. The ionization conditions probed by optical diagnostics using 
[ONIH versus [NII] /Ho and [SII]/Ha do not clearly indicate a (luminous) AGN, which at first glance my seemingly 
rule out the presence of a massive black hole. However, there are a number of potential reasons for this apparent 
discrepancy. First, [NII] and [SII] are not as sensitive to the hardness of the ionizing radiation field as [OI] and are 
therefore less reliable at identifying AGNs. Second, the massive black hole in Henize 2-10 is accreting at a tiny fraction 
of its Eddington luminosity (~ 107° L gaa) B8], which can the impact emission lines such that they do not look like 
those of luminous AGNs with higher Eddington ratios [43]. The apparent contradiction of emission line diagnostics can 
also arise when there is underlying star formation contributing to the spectrum. Finally, we note that there are other 
examples of radio-selected massive black holes in dwarf galaxies with enhanced [OI] that do not look like optical AGNs 
in the other diagnostic diagrams [35] [27]. 


We also find that the ionization conditions probed by strong emission line ratios at the location of the central source 
and along the filament are well described by theoretical models of shocks propagating in a high density medium 
(see Methods). Specifically, emission line ratios in the central ~50 pc along the EW slit orientation are consistent with 
models of a shock with a velocity of ~200 km s-1 traveling through a gas with an electron density of ~1000 cm~? and 
include a ‘precursor’ of ionizing photons that travel upstream from the shock front pre-ionizing the gas. The model 
shock velocities that match the strong emission line ratios agree with those obtained from direct measurements of the 
line widths from our kinematic study described above (e.g., [OII]), and the model gas densities are consistent with 
our measurements using the density sensitive line ratio [SII]6716/6731 (see Methods). Moreover, these conditions are 
well explained by an AGN-driven outflow mechanically exciting the interstellar medium in these regions in addition 
to the precursor component contributing to the photoionization of the interstellar medium. There is also evidence for 
shocked emission at the location of the eastern star-forming region, particularly from the secondary blue-shifted peak, 
in addition to photoionization from young hot stars (see Methods). Combining these results with our kinematic study 
indicates that the bipolar outflow generated by the central black hole is shocking the interstellar medium in the central 
regions of Henize 2-10, creating conditions that are favorable for positive AGN feedback [42] 4T]. 


AGN driven outflows have been discovered in a small sample of dwarf galaxies, though the discovery of a black hole 
outflow in Henize 2-10 provides the first example that is robustly spatially resolved. Moreover, Henize 2-10 differs 
from these other systems in several other ways. The majority of black hole driven outflows in dwarf galaxies have been 
found in galaxies with well-defined nuclei and optically-selected AGNs with relatively high accretion rates. These 
observations suggest the AGNSs play a role in heating and expelling gas in the galaxies and quenching star formation, 
a phenomenon known as negative feedback [33]. This is in stark contrast with Henize 2-10, which has an irregular 
central morphology, is intensely forming stars, and is experiencing positive feedback from a weakly accreting black 
hole that is luminous at radio, rather than optical, wavelengths. 


Indeed, “jet-mode” feedback is often associated with radio-loud AGNs accreting at low Eddington ratios, and attributed 
to the unbound nature of radiatively inefficient accretion flows [40]. Warm ionized gas outflows are observed to 
accompany the radio jets/outflows in some cases, particularly in young radio galaxies where nascent jets are expanding 
through the interstellar medium in the central regions of their hosts [40]. In these systems, the spatial extents of the 
warm outflows (traced by emission lines) are similar to the radio morphologies. A similar phenomenon is observed 
in Henize 2-10, which is most readily seen in a comparison between the radio emission detected with the Very Large 
Array and Paa emission detected by HST in the central few hundred pc of the galaxy [38]. However, an important 
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Figure 3: Visualization of the bipolar outflow model and star forming regions. Top left: Position-velocity diagram 
of the [OIT]5007 emission line along the EW slit where the Doppler shifts are relative to the systemic velocity of 
the galaxy (873 km s~ +) [21]. The data are well described by a simple precessing outflow model with a precession 
frequency of f ~ 5 revolutions Myr~! and a precession angle of 0 ~ 4° (see Methods). Bottom left: HST narrowband 
Ha-+continuum image showing that the EW slit position is aligned with the ionized filamentary structure connecting the 
massive black hole and site of recent star formation 70 pc to the east. White contours indicate compact radio emission 
from very long baseline interferometry [36], and show the location of the massive black hole. Bottom right: HST 
0.8 micron (~I-band) broadband image with the same field of view. Young star clusters in the eastern star-forming 
region are highlighted with an ellipse. Top right: Spectra of the eastern star-forming region (-70 pc) show strong broad 
emission lines with a secondary blue shifted peak most clearly visible in the [SII] emission. 


difference between Henize 2-10 and powerful young radio galaxies is that the emission lines are dominated by AGN 
photoionization in the more massive and luminous systems, whereas the extended emission line regions in Henize 2-10 
are dominated by star formation (enhanced/triggered by the black hole outflow). Therefore, Henize 2-10 may be a 
low-mass, low-power analog of young radio galaxies. 


1 Methods 


1.1 The Controversy — Massive Black Hole or Supernova Remnant: 


In recent years, evidence has been mounting for a massive black hole powering a low-luminosity active galactic nucleus 
(AGN) at the center of Henize 2-10 (38) (36) [37] [39], although a supernova remnant has been proposed as an alternative 
by some authors [16] [9]. As discussed in the main text, the radio and X-ray point source luminosities are consistent 
with both. A recent study argues for a supernova remnant based on their findings that the X-ray spectrum is better 
fit by a hot plasma model (typically used for supernova remnants) than a power-law model (typically used for luminous 
AGNs). However, the soft X-ray spectrum of the nuclear source in Henize 2-10 does resemble massive black holes 
accreting at very low Eddington fractions including Sagittarius A* in the Milky Way [2]. Another study using 
ground-based spectroscopy favored a supernova remnant origin for the central source based on a lack of any AGN 
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Observations massive supe = Notes and References 
black hole? remnant? 
Nuclear, non-thermal + Nuclear, non-thermal point source detected with the VLA (Johnson et al. 2003; Reines et al. 2011) 
radio point source with yes yes + VLBI observations confirm non-thermal origin and constrain the nuclear radio source to be <3 x 1 
la~4x 10% erg/s pe. If SNR, size and luminosity imply t must be very young (Reines & Defer 2012). 
Spatially coincident + Nuclear X-ray point source detected with Chandra that is spatially coincident with radio source. If 
nuclear X-ray point massive BH, X-ray luminosity implies highly sub-Eddington accretion, which is consistent with soft 
h yes yes X-ray spectrum. An X-ray binary is ruled oul by the high luminosity of the compact radio source 
source with Lx ~ 10% (Reines ot al, 2016) 
erg/s + X-ray spectrum also consistent wath SNR (Hebdar et al. 2019) 


X-ray light curve from 160 ks Chandra observation exhibits moderately sign#icant variability (a 
X-ray variability on hr- factor of ~2: Reines et al. 2016). Probability of non-variable Sight curve producing the signal is 3% 
long timescales yes no (Hebbar et al. 2019) 

+ Tentative detection of ~9-hr periodicity in the variable light curve (Reines et al. 2016) 


+ Stellar velocity dispersion map from near-IR IFU observations reveals enhancement of ~7 km/s at 


Enhanced stellar position of nuclear source, consistent with a BH mass of ~1.5 x 10® Mu» (Riffel et al. 2020). This is 
velocity dispersion at in Ene with expectations given the stellar mass of Henize 2-10 (~10"° Mun, Nguyen et al. 2014) 
location of nuclear yes no + No detectable star cluster at location of the central source (Reines et al, 2011). This is in contrast to 


a weaker, more extended off-nuciear VLBI source, which is co-spatial with the most luminous super 


aes star cluster in Henize 2-10, and is a tkely SNR (Reines & Deller 2012). 


« Enhanced near-IR [Fe il/Gry ratio consistent wiih accreting BH of supernova driven shocks (Rite! 
@t al. 2020; Crosci et al. 2010) 
Gas excitation and 9 + Near-iF Hz line ratios consistent with AGN photoionization (Rittel et al. 2020) 
ionization conditions yes . + Optical line ratios from ground-based spectra dominated by star formation (Cresci et al. 2017) 
HST/STIS observations show non-stelar emission-line ratios involving [Ol]. SNRs have log [Silja 
> -0.5, which is not seen in the STIS spectrum of the nuclear source (this work) 


* Spectroscopy at high spatial resolution shows broadened near-IR and optical emission lines at the 
location of the central source (Riffel et al. 2020 and this work). in particular, (06300 has a FWHM 
of ~500 km/s indicating an outflow from the central source (this work), Very young SNRs typically 

Gas kinematics and have broad lines with widths on the order of thousands of km/s 

morphology yes no * HST/STIS spectroscopy along a ~150 pc-long ionized gas filament centered on the nuclear source 
(seen in Ha imaging with HST) exhibits a coherent sinusoid-like pattern in the Doppler shifted 
velocities, The position-velocity data is well-tit by a simple mode! of a processing bipolar outflow 
from a massive BH, and is incompatible with a SNR (this work) 


Figure 4: (Extended Data Table 1) Summary of observational results regarding the nature of the nucleus in 
Henize 2-10. A highly sub-Eddington massive black hole is consistent with all the available observations including 
new results presented here, while a supernova remnant is not. 


ionization signatures [9]. However, the ground-based observations used in that work had an angular resolution of 
~0.7", which is not sufficient to cleanly isolate the weakly accreting black hole from nearby young (<5 Myr) massive 
(M, ~ 105 Mo) star clusters that dominate the line ratios at this relatively course angular resolution. There are other 
observational results to consider regarding the origin of the nuclear radio/X-ray source in Henize 2-10. For example, 
a recent study using adaptive optics integral field spectroscopy provides evidence for gas excited by an AGN and 
an enhanced stellar velocity dispersion at the location of the nuclear source consistent with a ~ 10° Mo black hole, 
favoring the low-luminosity AGN interpretation [39]. There is also evidence for moderately significant variability on 
hour-long timescales in the X-ray light curve, which is incompatible with a supernova remnant [37] [16]. Moreover, it is 
reasonable to expect that Henize 2-10 hosts a massive black hole since its overall structure resembles an early-type 
galaxy (albeit with a central starburst) and its stellar mass may be as high as M, ~ 10!°Mo [BO], a regime where the 
black hole occupation fraction is near unity [15]. The central starburst complicates the identification of the weakly 
accreting black hole, yet a variety of multiwavelength observational results taken collectively strongly support its 
presence. These results are summarized in Figure[4](Extended Data Table 1). A highly sub-Eddington massive black 
hole is consistent with all of the observations, including the new work presented here, while a supernova remnant is not. 
The present study not only adds to the evidence for a massive black hole in Henize 2-10, it also demonstrates that a 
bipolar outflow from the black hole is enhancing/triggering star formation in its vicinity. 


1.2 STIS Observations and Data Reduction: 


Spatially resolved spectroscopic observations of the nuclear regions of the dwarf starburst galaxy Henize 2-10 were 
obtained using the Space Telescope Imaging Spectrograph (STIS) instrument on the Hubble Space Telescope (HST). 
We obtained observations with two slit orientations. The first is aligned with the quasi-linear ionized gas structure 
identified by Reines et al.[38] and covers the central radio/X-ray source and the bright knot of ionized gas to the east. 
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Figure 5: (Extended Data Figure 1) Raw 2D spectra showing the [O1]6300 emission line at the location of the 
nucleus in the EW slit orientation. The location of the nucleus is indicated by white circles and the two images 
correspond to the two dithered sub-exposures. 


We refer to this as the East-West (EW) orientation. The second slit orientation was placed perpendicular to the EW 
observation at the location of the central radio/X-ray source. We refer to this as the North-South (NS) orientation. The 
candidate AGN itself was too faint to acquire directly, therefore we used a target acquisition with an offset from a bright 
point source 7.9" to the southeast. Spectra were taken with the G750M and G430M gratings providing medium spectral 
resolution (R ~ 5000-6000) coverage of key optical emission lines. The central wavelengths were set at 6581 A and 
4961 A for the G750M and G430M gratings, respectively. At each slit orientation, two orbits were spent in G430M and 
one orbit in G750M. The observations were taken with a two-point dither pattern with CR-SPLITS (multiple exposures 
taken to aid in cosmic ray rejection) at each position to help eliminate cosmic rays. The calibrated dithered images were 
combined and have a spatial resolution of ~0.1", which corresponds to a physical scale of ~4 pc at the distance of 
Henize 2-10 (~9 Mpc). 


1.3 Emission Line Fitting: 


Before we fit emission lines in the spectra, we first modeled and removed the continuum. The reduced spectra were 
continuum-subtracted by masking emission line regions and fitting a low order polynomial to the continuum in each 
row in the spatial dimension. A low order polynomial was used given the lack of absorption features in the spectra. 
We did, however, consider the potential impact of stellar absorption lines on our measurements and found that the 
absorption line strengths are negligible compared to the emission line strengths. Scaling a Starburst99 model for a 
4 Myr stellar population (see Stellar Ages section below) to our observed spectra, we find that the flux of Ha absorption 
is smaller by a factor of 61 than the Ha emission and the flux of the H£ absorption is smaller by a factor of 7 than 
the H8 emission at the location of the central source. Accounting for this absorption has a negligible effect on the 
line ratios of the nuclear source and does not impact the classifications based on the diagnostic diagrams. Once the 
spectra were continuum-subtracted, we fit each emission line of interest with a linear combination of Gaussian profiles 
to characterize the flux and estimate kinematic properties along the spatial dimension of each slit. The fitting was 
done using Imfit [29], a non-linear least squares curve fitting package in Python. We fit each emission line with up 
to two Gaussian components when needed. To determine if a second Gaussian component is warranted, we require 
that the flux of both components be greater than the 3a error of the flux. This process is performed row by row in 
the spatial direction along each slit for emission lines of interest. During this process we fit the HG, [OII]5007, Ha, 
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Figure 6: (Extended Data Figure 2) Combined 2D spectra showing the [01]6300 emission line at the location 
of the nucleus in the EW slit orientation. Same as Figure|5|but showing the reduced 2D image with the dithered 
sub-exposures combined. 


[NII]6548/6583 and [SII]6716/6731 emission lines. We fit the Ha and [NII] lines simultaneously, fixing the spacing 
between [NII] Gaussian components to their corresponding component in the Ha emission line to laboratory values. 
Additionally, we tie the widths of [NII] components and fix the flux ratio of the [NII] lines to the laboratory value of 
1:2.96. Similarly, the two [SII] emission lines are fit simultaneously with the spacing between Gaussian components of 
the two lines held fixed and the widths of the Gaussian components are tied together. We also fit [OI]6300 in the spectra 
of the nuclear source and the eastern star-forming region, but the line is too weak to be detected all along the slits. 
Since the [OI]6300 line has a complex profile at the location of the black hole, with possible double peaked narrow 
lines and a much broader component than the other emission lines, we confirmed that this was not due to an artifact 
in the data. In Figure|5| we show close-up views of the raw 2D spectra along the EW slit position. The two images 
correspond to the two dithered sub-exposures offset by 7.5 pixels. The broad [OI] line at the location of the central 
source is seen in both images at different positions on the detector (indicated by white circles). Note that the locations 
of hot pixels do not change between the two images. In Figure|6|we also show the final reduced 2D image with the 
sub-exposures combined. The broad, double peaked nature of the [OI] line is clearly visible. We also note that [OI] 
is similarly broadened in the nuclear spectrum extracted from the NS slit position, although there is not an obvious 
double-peaked narrow line component (see Figure 2 in the main paper). In any case, broadened [OI] is clearly detected 
at the location of the nuclear source in both slit positions indicating an outflow on the order of 500 km s~!. 


1.4 Gas Density: 


We estimate the electron density, ne, using the density sensitive line ratio [SII]6716/6731 [32]. This ratio is sensitive 
to electron densities in the range of ~10?-104 cm~’. Along the EW slit orientation, we find a range of ne ~ 
1025 — 10*cm~3, indicating a relatively high-density gas (see Figure|7). These density estimates are in general 
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Figure 7: (Extended Data Figure 3) The electron density, ne, along the EW slit orientation. We measure the 
electron density along the EW slit from the ratio of [SII]6716/[S1I]6731 and find the electron density ranges from 
~ 1025 — 104cm~? , which is within the range the [SII] ratio is sensitive to density. The high densities are consistent 
with those predicted by optical emission line diagnostics derived from the Allen et al. shock models. 


agreement with the gas densities predicted by the Allen et al. shock/shock+precursor models in the central regions 
of Henize 2-10 as described in the next section. 


1.5 Emissin Line Diagnostics - Photoionization and Shock Models: 


To understand the ionization mechanisms in the central regions of Henize 2-10 we compare our emission line mea- 
surements in various regions to photoionization and shock models. In addition to the central radio/X-ray source, we 
identified 7 regions of interest that are shown in Figure [8] and serve to provide a spatially resolved picture of the 
kinematics and ionization conditions in the central regions of Henize 2-10. The extraction regions taken along the EW 
slit orientation were chosen to correspond with emission features seen in the Ha and I-band imaging from HST (young 
star clusters, knots of ionized gas) as well as features seen in the STIS spectroscopy (broad emission, double peaks). 


We first utilize the standard emission line diagnostic diagrams described by Baldwin et al. and Veilleux et al. 
that have been expanded upon and summarized in Kewley et al. [19]. An accreting BH will produce a much harder 
continuum than is emitted by hot stars, and these diagrams take advantage of this fact by comparing strong emission 
line ratios that are close together in frequency to mitigate reddening effects. In this study we employ widely used 
emission line diagnostic diagrams that take [OII]/H( versus [NII]/Ha, [SII]/Ha, and [OI]/Ha (see Figure p}. In the 
[NII]/Ha diagram, line-emitting galaxies separate into a V-shape with star forming galaxies occupying the left most 
plume while AGNs occupy the right branch of galaxies. These regions are quantified by an empirical division between 
HII regions and emission from AGNs developed by Kauffmann et al. [17]. The “composite” region between this 
empirical division and the theoretical maximum starburst line from Kewley et al. indicates there is likely significant 
emission form both HII regions and AGNs. Like the [NII]/Ha diagram, the [SII]/Ha and [OI]/Ha diagnostics provide 
diagnostics for differentiating between emission from HII regions and AGNs. These two diagrams add a dividing line to 
distinguish between emission from Seyferts and Low Ionization Nuclear Emission Regions (LINERs). LINER emission 
can be generated both by shocks and very hard AGN spectra and determining the primary ionization mechanism can be 
complicated [B]. It should be noted that while these diagnostic diagrams are useful for identifying regions dominated by 
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Figure 8: (Extended Data Figure 4) The spatial extraction regions taken along the EW slit orientation. We place 
these regions on optical emission line diagnostic diagrams (Figures|9}11). Top panel: the extraction regions are shown 
on the narrow band Ha + continuum image from HST to highlight the ionized gas features that several of the spatial 
extractions probe. Bottom panel: the extraction regions are shown on the archival 0.8 micron HST image, showing 
young star clusters that the EW slit orientation passes through. 
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Figure 9: (Extended Data Figure 5) Narrow emission line diagnostic diagrams showing various extraction re- 
gions along the EW slit orientation (see Figure [8). The nucleus (yellow point) falls in the Seyfert region of the 
[OI]/Ha diagram. The young star-forming region ~70 pc to the east of the low-luminosity AGN is depicted with a blue 
triangle and star for the primary emission line component and the blue-shifted secondary component, respectively. [OI] 
is not detected in all of the regions. 


luminous AGNSs, they have limitations and can yield ambiguous results for (or completely miss) massive black holes 
accreting at very low Eddington ratios such as the one in Henize 2-10. Indeed, non-stellar ionization is clearly indicated 
in the [OI]/Ha diagram at the location of the nuclear source, but not so for the other diagnostic diagrams (see discussion 
in the main paper). Figure|9]shows the [OIII]/H( versus [NII]/Ha, [SII]/Ha, and [OI]/Ha diagnostic diagrams for the 
various extraction regions along the EW slit. 


In addition to the diagnostic diagrams discussed above, we also investigate whether the ionization conditions seen in 
the central regions of Henize 2-10 can be explained by mechanical excitation from shocks. To investigate this, we 
employ ionization models of shock and shock+precursor emission developed by Allen et al. [I]. These models provide 
emission line fluxes for ionization from a pure shock (possibly driven by an outflow from an AGN or by regions of 
intense star formation), where the gas is collisionally ionized, or a shock+precursor where ionizing photons produced in 
the shock-heated gas travel upstream and ionize the gas before the shock reaches it. We explore models with a variety 
of electron densities (0.01-1000 cm~), shock velocities (100-600 km s~!) and transverse magnetic field strengths 
(0.01-32 uG). These are shown in Figures[I0]and[I1] 


The emission line ratios from the nuclear source are best described by the shock+precursor models with a low shock 
velocity (100-250 km s71), a high-density gas (n = 1000 cm~3), and a low transverse magnetic field parameter (b = 
0.01 — 1 uG) (Figures[I0]and[I 1). Shock+precursor models are thought to be a good description for AGN+outflow 
emission. Along the filament (extraction regions 2-5), the line ratios are explained well by a low velocity shock or 
shock+precursor model (~200 km s~') in a high density (ne ~ 1000 cm~?) gas with a transverse magnetic field 
parameter in the range of 1-10 uG. This is consistent with a scenario where the central black hole is driving a bipolar 
outflow that shocks the gas and dominates the ionization conditions along the filament. 


At the region of intense star formation located ~70 pc to the east of the black hole, we observe strong emission lines 
including a secondary, blue-shifted kinematic component. To properly fit the emission lines in this region, an additional 
Gaussian component is added (see Figure ]in the main paper). The separation of this component is determined using 
the [SII] emission line region where the secondary peak is most clearly resolved. We find that the secondary peak is 
offset from the primary peak by 154 km s~+, and we use this value when fitting the other emission line regions where 
the secondary peak is not as well resolved. When comparing to shock(+ precursor) models, we find that the [NII]/Ha 
diagram indicates a low shock velocity (100-250 km s71) in a high-density gas (n = 1000 cm~?) with a low transverse 
magnetic field parameter (b = 0.01 — 1 uG) (see Figures [I0]and [1 1}. The [SII]/Ha diagram shows the primary emission 
peak is inconsistent with emission from shocks or shocks+precursor models. This indicates that the primary emission 
peak at this location is primarily due to star formation. The secondary emission peak is consistent with shock+precursor 
models for the low velocity, high density conditions, indicating that this kinematically distinct emission component is 
dominated by a shock+precursor from the AGN-driven outflow. 


Emission from extraction regions 6 and 7 (i.e., the western star-forming region) is not consistent with any shock or 
shock+precursor models, which is in agreement with their location in the HII region of the Baldwin—Phillips—Terlevich 
(BPT) diagram. The line ratios are dominated by star formation in this region. 
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Figure 10: (Extended Data Figure 6) Optical emission line diagnostics from the shock and shock+precursor 
models with varying gas density. We place the spatial extractions from the EW slit orientation shown in Figure[8]on a 
grid of shock excitation models (presented in Allen et al. [T]) with varying gas density (n = 0.01-1000 cm~) and shock 
velocity (v = 100-600 km/s). We fix the transverse magnetic field to be b = 1 uG and the assume solar metallicity. 
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Figure 11: (Extended Data Figure 7) Optical emission line diagnostics from the shock and shock+precursor 
models with varying magnetic field. The models (presented in Allen et al. [T]) are shown as a grid with dashed blue 
lines indicating constant shock velocity and dashed black lines indicating constant transverse magnetic field. For these 


models, the density is fixed to n = 1000 cm~? and the transverse magnetic field parameter is allowed to vary from b = 
0.01-32 uG. 
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1.6 Star Cluster Ages: 


We estimate the ages of the young stellar clusters that fall within the EW slit from their Ha and H£ equivalent widths. 
To ascertain ages from these equivalent widths we employ simple stellar population (SSP) models from Starburst99 
22]. We use models from Version 7.1 with solar metallicity (appropriate for the central regions of Henize 2-10 [25]), 
an instantaneous burst of 104 Mọ with a Kroupa IMF (0.1 — 100 Mo), the Geneva evolutionary tracks with high mass 
loss and the Pauldrach-Hillier atmospheres. 


At the location of the young stellar clusters in the eastern star-forming region (region 1 in Figure[8) we find an equivalent 
width of 478 A and 70 A for Ha and H8 respectively. These both give stellar population age estimates of ~4.3 Myr, 
which is in good agreement with previous estimates of the ages of other young star clusters in the region [6]. The ages 
of these clusters are larger than the crossing time for the AGN-driven outflow (~0.3 Myr), based on the minimum 
outflow velocity measured from emission line spectra (~200 km s71) and the distance between the AGN and the eastern 
star-forming knot (~70 pc). Therefore, the timescales allow for the AGN-driven outflow to have triggered/enhanced the 
formation of star clusters in the Eastern star-forming knot. 


The EW slit orientation also passes through a young star cluster in region 3. At this location we find equivalent widths 
of 212 A and 41 A for Ha and H8 respectively, both indicating an age of ~5.2 Myr for the stellar cluster. Finally, 
region 6 in the western star-forming region has Ha and H£ equivalent widths of 1092 A and 196 A, respectively. These 
equivalent widths indicate the stellar clusters have ages <3 Myr. 


1.7 Bipolar Outflow Model: 


Here we provide a derivation of the model used to describe a precessing bipolar outflow emanating from the central 
radio/X-ray source, which can explain the coherent velocity structure seen in the central ~120 pc of the EW orientation 
observations. In this model we align the EW slit orientation with the z-axis and assume the outflow precesses about this 
axis with a small angle 0 and an angular precession frequency w. If the gas being ejected by the outflow has velocity vo, 
and we orient the x-axis to be in the direction of the observer, then the radial (Doppler shifted) velocity seen at the 
location of the AGN as a function of time will be given by 


vu,(t) = vz(t) = vo sin (0) sin (wt + 7), (1) 


where y represents a phase shift that accounts for small asymmetries in the outflow profile (see Figure [12] for an 
illustration of our model). 


In order to find the radial velocity as a function of distance (z) along the slit axis, we must consider what angle the 
outflow made with the (line-of-sight) x-axis when the gas at distance z was emitted. Since this angle is time dependent 
as the outflow precesses, the line-of-sight velocity of the gas will depend on the orientation of the outflow at some time 
to in the past. The time that has passed since the gas at distance z was ejected by the outflow is determined by the z 
velocity of the gas. Due to the symmetry of the model about the z-axis, the z component of the gas velocity will be 
time independent and only depend on the angle of the outflow with the z-axis, 


Uz = vo cos (8). (2) 
The time, tg, for gas to reach a distance z along the slit is then given by 


zZ 
o= acos) a 


We are then able to find the an expression for v, (z) by evaluating the expression for v, (t) at the time —to: 


v,(z) = vo sin (0) sin (. — 7°) (4) 


vo cos (0) 
To fit this model to the data we require a rough estimate of the bulk outflow velocity, vo. We estimate this parameter 
using W80, the velocity interval containing 80% of the line flux, of the broad emission seen at the location of the 
candidate AGN. We find W80 ~ 200 — 500 km s~! based on measurements of the [OIII]5007 and [O1]6300 lines at 


the location of the candidate AGN. This allows us to determine the best-fit angle of precession, 0, and the frequency of 
precession, f = w/2r, to be 
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0 = 2.4° — 6.1° (5) 
f = 3.0 — 7.5 revolutions Myr’ (6) 


where the larger angle of precession and smaller frequency of precession corresponds to lower outflow velocities. We 
find consistent results when using the Doppler shift profile of Ha, H8 and [OII] emission lines to fit the model derived 
above (the results using the [OII] emission line is shown in Figure Bp. The Doppler shift profile can be coherently 
traced out to 50-60 pc on either side of the candidate AGN, most definitively out to the bright eastern star forming 
region after which the Doppler shifts of the emission lines are influenced by the bright young stars and then shown 
no coherent pattern further along the slit. The coherent velocity structure seen on the scale of 100 pc is not consistent 
with young supernova remnant as the compact radio/x-ray source in the central regions of Henize 2-10, which provides 
further motivation that a low luminosity AGN is driving the outflow. 


These results are roughly consistent with jet parameters derived in other studies where precessing or reorienting jet 
models have been applied. The long precession period we observe (~200,000 years) is shorter by a factor of a few than 
those seen predicted by Dunn et al. [10], Nawaz et al. and Cielo et al. [[7] but longer by a factor of a few than those 
predicted by Gower et al. when jet precession is invoked to explain the complex bending and knotting seen in large 
radio jets. 
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V 


outflow 


Figure 12: (Extended Data Figure 8) A diagram of the toy model of the bipolar outflow generated by the low- 
luminosity AGN in Henize 2-10. Our simple model depends on the outflow velocity of the ionized gas (Voutflow), the 
angle the outflow makes with its precession axis (0) and the angular frequency with which the outflow precesses (w). 
Similar models have been used to describe the bending seen in large radio jets [14] [10]. 
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